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Abstract 

We present the results of molecular dynamics modeling on the structural properties 
of grain boundaries (GB) in thin polycrystalline films. The transition from crystalline 
boundaries with low mismatch angle to amorphous boundaries is investigated. It is 
shown that the structures of the GBs satisfy a thermodynamical criterion suggested 
in ||. The potential energy of silicon atoms is closely related with a geometrical 
quantity - tetragonality of their coordination with their nearest neighbors. A crossover 
of the length of localization is observed. To analyze the crossover of the length of 
localization of the single-electron states and properties of conductance of the thin 
polycrystalline film at low temperature, we use a two-dimensional Anderson localization 
model, with the random one-site electron charging energy for a single grain (dot), 
random non-diagonal matrix elements, and random number of connections between 
the neighboring grains. The results on the crossover behavior of localization length of 
the single-electron states and characteristic properties of conductance are presented in 
the region of parameters where the transition from an insulator to a conductor regimes 
takes place. 
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Recently it was demonstrated that thin polycrystalline silicon films are promising 
materials for future room temperature single-electron devices |T|-H • Main reasons 
which make this material so attractive are the following: I). Usually, the film's thickness 
varies from 1 to 5 nm, and the average lateral grain size is 10 nm or less. In this case, 
the energy of an electron in a single grain, E e , is bigger than the thermal energy 
even at room temperature, E e > T = 300° K. II). It is believed that in these films, 
the characteristic resistance, i4, of the potential barriers between the grains is big 
enough, > R c = h/e 2 ~ 25kQ. If both these conditions are satisfied, an electron is 
strongly localized in the grain. At the same time, one can regulate (up to some extent) 
the electron conductivity in these films by varying the gate voltage, and creating a 
current channel. Those electrons which are stored in the grains (storage dots) create 
a Coulomb repulsion for those electrons which are involved in the current channel. 
This allows one to implement memory operations in these films at room temperature 
using a Coulomb blockade effect. Different implementations of these ideas have been 
discussed, for example, in To satisfy both above mentioned conditions, the 

boundaries between the nano crystalline grains in these films play an important role. 
For example, one of the most important characteristics of the electron transport in 
polysilicon films is connected with the distribution of crystalline and amorphous grain 
boundaries (GB) ||. At room temperature, the main factor which determines the 
structure of a GB is a mutual mis-orientation of the neighboring crystalline grains. At 
present there does not exist a consistent theoretical approach for the description of the 
GBs in thin polycrystalline films. 

I. Molecular dynamics (MP) simulation of grain boundaries at room temperature. 

The initial configuration of the polysilicon film was generated using the Voronoi 
triangulation algorithm: first, the centers of the grains were defined; then each grain 
was obtained by filling the space, nearest to the center of the grain, with the diamond 
lattice at a chosen orientation. The atoms with very high potential energies were 
eliminated from the GBs. The film had the dimensions 16.7 x 16.7 x 3.1 nm, and 
contained 44174 atoms in 12 grains (Fig. 1). The majority of the grains in our model 
were oriented randomly, so they had both tilt and twist mis-orientation. Four grains, 
namely the 3,4,7 and 8-th were generated so that the mis-orientation tilt angle between 
the 3-rd and 4-th grains was a 3 _ 4 = 11.0°; between 4-th and 7-th grains a 4 __ 7 = 11.0°; 
and between 7-th and 8-th grains, a-j_% = 15.0°. The GB 3 _^ and the GB 7 _ S are close 
to (l,l,l)-interface, and the GB^_ 7 is close to (0,0,l)-interface. The (l,l,l)-planes 
in these four grains are perpendicular to the plane of the figure. The average linear 
size of a grain in our model, d gr ~ 4.8 nm, is of the same order as for the grains in 
the experimental polycrystalline silicon films obtained recently by thermal annealing 
of amorphous Si deposited by decomposing silane on SiC>2 substrate ||]. Periodic 
boundary conditions were used in the X and Y directions, and free boundary conditions 
were used in the Z - direction. It is shown that the potential energy of the atoms 
closely correlates with the degree of their tetragonal coordination with the nearest 
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neighbors. To estimate the quality of the tetragonal coordination of the j-ih atom, we 
calculated the value, Tj = J2i<k(lji ~ (j'fc) 2 )/!^, which is called the tetragonality 0. 
Here Ijk is the length of the k-th edge of the tetrahedron of general shape formed by 
the four nearest neighbors of the j-th atom; lj is the average length of the edge. By 
definition, Tj is equal to zero for an ideal tetrahedron (as it is the case of ideal diamond 
structure), and increases if the shape of the tetrahedron is distorted. The value, Tj, is 
correlated with the dispersion of the edge lengths of the tetrahedron. According to the 
criterion suggested in ||, the spatial figure formed by four points can be recognized as 
having the "good tetrahedral shape" if Tj is less than or equal to T^°> = 0.018. The 
value of tetragonality of an atom is very sensitive to the number of atoms in the first 
coordination shell, z. In the bulk amorphous phase of silicon, for which z a ~ 4.04, the 
average tetragonality is equal to T am = 0.015 which is slightly less than the critical 
value T^ c \ In liquid silicon (in which the first coordination shell consists in average 
of zi q ~ 4.46 neighbors), the average value of tetragonality is, < T >^~= 0.034. In 
Table I, we present the data for the potential energy per atom, density and average 
tetragonality, for the longest GBs (for which it was possible to make measurements 
with reasonable accuracy). Densities of the tilt crystalline GBs between the 3, 4, 7 
and 8-th grains are practically the same, or slightly smaller, than the density of the 
bulk crystal, p cr . The densities of other GBs are systematically larger than p cr . In 
covalent materials with friable structure, such as silicon and germanium, some GBs, 
indeed, contract. This was observed in X-ray diffraction experiments for (1,0,0) twist 
GBs in Ge M and in MD simulation of the (l,l,l)-twist boundaries in silicon [pUfl . It 
is remarkable that the contraction is observed for the twist boundaries. In such GBs, 
the covalent bonds are most probably stretched rather than squeezed, and the system 
demonstrates a tendency to restore their lengths at the cost of a reduction of the GB's 
volume. MD simulations show that in our model of polysilicon film, the majority 
of the GBs between the grains with random mutual orientation are disordered. The 
grains with small mutual mis-orientation are connected by crystalline GBs. The GBs 
with "medium-angle" mis-orientation have a complex inhomogeneous structure. These 
boundaries consist of crystalline connections interspersed with disordered regions. In 
disordered regions, the misfit between the crystal structures of the neighboring grains is 
compensated by non- crystalline arrangements of atoms, such as 5- and 7- fold rings. In 
general, the structure of the GBs satisfies the thermodynamical criterion suggested in 
|J. The potential energy of the silicon atoms closely correlates with the tetragonality, 
T, of their coordination with their nearest neighbors. The majority of atoms in the 
polysilicon grain boundary are well coordinated tetrahedrally, even if they are arranged 
in a non-crystalline manner in the high-angle GB. From the observed small values 
of the average tetragonality, it follows that the high-angle disordered GBs can be 
characterized as amorphous. 

The results obtained are important for better understanding of structural properties 
of thin polycrystalline silicon films, which have been used recently for memory devices. 
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However, for obtaining more reliable data, larger polysilicon simulations should be 
considered. 

II. Crossover behavior of the localization length of single-electron states, 
and the properties of conductivity. 

To model the properties of conductive electrons in thin polycrystalline silicon film at 
low temperature, we use an approach based on the two-dimensional tight-binding An- 
derson model. Each nanosilicon grain (with the coordinates (m, n) in (X, Indirections) 
is considered as an individual dot. The characteristic energy of quantization in the Z- 
direction, E z , for the lowest level, is of the order, E z ~ 1 — 300 meV [|I|, E|. We 
assume that the total one-site electron energy is E m ^ n =< E > +5E m ^ ni where 5E m n 
is a random variable (as the thickness of the film fluctuates). At low temperature, 
the electron conduction is governed by quantum tunneling between the dots, which 
we describe by the matrix elements, Vj^' n \ Below we consider only the simplest 
case of the neighboring transitions: V^™ ±1,n±1 ) =< V > +5V mjn , where < V > is 
the average value of the off-diagonal matrix elements, and SV m>n describes the fluctu- 
ations related with the random characteristics of the potential barriers between the 
neighboring dots. Let us now estimate the number of conduction electrons which oc- 
cupy an individual dot. The maximum value of a 2D electron density in the film can 
be chosen as: n e ~ 10 11 — 10 12 cm -2 . Consider a subsystem of conduction electrons 
in a polysilicon film as a two-dimensional one, with the size of the active region 100 
nmxlOO nm |T|, Qj. Then, the number of conduction electrons, Af, in this region is of 
the order M = 10 — 100. If the average size of the grain (dot) in (X, Y) plane is of the 
order lOnm, then one has, on average, 100 dots in the active region. So, the average 
number of conductive electrons in a dot is of the order, 0.1-1. In this paper, we con- 
sider the conduction electrons as noninteracting. To describe the electron subsystem 
at low temperature, we use a 2D Anderson tight-binding model with the Hamiltonian 
TC — J2i Ei\l >< l \ + J2i k Vi,kV X k\i where E\ describes the diagonal disorder. Fig. 2 
shows the dependences of the average length of localization, < L >, and the dispersion 
of the length of localization, D(L) =< (L— < L >) 2 >, of the single-electron states, 
\[>W(m, n), as a function of the average value of the matrix element < V >. (The 
size of the lattice is (M, AT) = (18, 12). Zero boundary conditions were chosen.) In Fig. 
2, < E >= 50 meV; £ m , n G [50 - 12.5; 50 + 12.5] meV; W m , n G [0.9; 1.1] < V > 
meV. For each state, % = 1,...,MN, the length of localization was introduced as 

L» = (L« + L«)/2. We calculated L® we as L« = sJyS^Li (m - ? I * ijn, n) \ 2 , 

where = J2m,n=i m \^(. m i n )\ 2 - (Similar expressions were used to calculate L$ .) 
As one can see from Fig. 2, the fluctuations of the average localization length, D(L), 
exhibit a characteristic maximum at < V >~< V c >= 3.5 meV. For < V ><< V c >, 
the system is essentially insulating. For < V »< V c >, the system exhibits metallic 
properties. We are presently investigating this crossover behavior in connection with 
the metal-insulator transition (MIT). Finally, in Fig. 3, we present the characteris- 
tic results of numerical simulations of the dimensionless conductance in this system, 
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<y = Yj \T a ,j3\ 2 -, where T a ^ is the transmission amplitude from channel a to channel f3. 
The method we used is based on a Green's functions approach flr2"|. The lattice size 



is (11,11), with 11 channels in the left and the right electrodes. In Fig. 3, the hori- 
zontal axis is the dimensionless value of the diagonal disorder, 5E. The dimensionless 
off-diagonal matrix element was chosen equal to V = 1 in the active region and in the 
electrodes. The results shown in Fig. 3 correspond to the dimensionless Fermi energy 
Ep = 0. In Fig. 3 the data indicated by o (a) corresponds to the randomization of only 
two boundaries, (n=l,m = l,..., M) and (n = N, m = 1, M). The data indicated 
by + (b) corresponds to the randomization of the whole active region. One can see that 
increasing the diagonal disorder leads to a crossover behavior of the conductance in 
the case when only the boundaries are randomized. The results presented in this paper 
can be useful for performance the memory devices based on thin polycrystalline silicon 
films. Our further work is related to studying the influence of the interfaces Si/Si02 
on the conductivity in polycrystalline silicon films, and on the noise component. 
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U (eV/ at), ±0.04 


p (g/ cm 3 ), ±0.02 


< T >, ±0.0003 


GB 3 _ 4 


-4.24 


2.29 


0.005 




-4.18 


2.31 


0.014 


GB 7 _ 8 


-4.20 


2.34 


0.011 


GB 3 _q 


-4.19 


2.39 


0.010 


GBq- 7 


-4.16 


2.39 


0.013 


GB±_i 


-4.17 


2.41 


0.013 


G-Bi_5 


-4.18 


2.38 


0.011 


GB 5 - 2 


-4.09 


2.36 


0.017 


GBg_i 


-4.21 


2.37 


0.012 


C-Bio-n 


-4.15 


2.38 


0.015 




-4.14 


2.38 


0.017 


a — Si 


-4.15* 


2.29* 


0.015* 


cr — Si 


-4.335 


2.324 






Table I: Potential energy, density and average tetragonality of grain boundaries 
in the polysilicon model. *Data for a-Si obtained from the model |7J 
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Figure Captions 



Fig. 1. The (X, y)-projection of the polycrystalline film. 

Fig. 2. Demonstration of the crossover behavior of the average localization length 
of the single-electron states. Dependence of the average localization length < L > 
(a), and of the dispersion D(L) (b) on the average value of the matrix element, 
< V >. (M, N) = (18, 12); < E >= 50 meV; SE = 25 meV; V =E [0.9; 1.1] < V > 
meV. 

Fig. 3. Dependence of the conductance, Z)T aj( g, on the disorder parameter SE. 
The lattice size (M, N) = (11, 11). 11 channels in the left and right electrodes were 
used. Dimensionless parameters: V = 1; Ep = 0; (a) Two of the boundaries are 
randomized: (n — 1, m — 1, ...,M), and (n — N, m — 1,...,M); (b) The whole 
active region is randomized. 
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Fig. 1 
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